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Abstract 

We calculate the evolution of the prompt intrinsic Poynting power generated by merging black 
holes. Orbiting black holes induce rotation of the space-time. In a presence of magnetic field 
supported by an accretion disk outside of the orbit, this results in a generation of an electromagnetic 
outflow via the Blandford-Znajek-type process with luminosity Lem ~ G 3 M 3 B 2 / (c 5 R or b) and 
reaching a fairly low maximum values of Lem = 10 37 — 10 39 m6 erg s -1 (rriQ is the masses of 
black holes in millions of Solar mass) at the time of the merger. Dissipation of the wind power 
may produce two types of observed signatures: a highly variable collimated emission coming from 
the internal dissipation within the jets and a broad-band near-isotropic emission generated at the 
termination shocks. 



PACS numbers: 04.30.Tv, 95.85. Sz 



I. INTRODUCTION 

Direct detection of gravitational waves of the associated merging events of compact ob- 
jects is one of the prime goals of modern astrophysics. Observations of the corresponding 
electromagnetic signal is most desirable, as it will provide crucial information on the loca- 
tion and the physical properties of the event. In case of the merger of stellar-mass compact 
objects (neutron stars and/or black holes), the short Gamma Ray Bursts may be the cor- 
responding electromagnetic event pp. On the other hand, the detection of the merger of 
supermassive black holes in the centers of galaxies by the Laser Interferometer Space An- 
tenna, is an exciting possibility [2, 3J. There are clear observational evidence for recent 
galaxy mergers jl] and direct observations evidence of binary black holes [5] . 

So far, the studies of a possible prompt electromagnetic signatures of the black hole 
mergers were attributed mostly to the perturbations that the black holes induce in the sur- 
rounding gas, [6H9]. The resulting electromagnetic signal is then subject to great uncertainty 
and naturally depends on the complicated non-linear fluid behavior of the system. One of 
the problems is that there should be little gas inside the orbit of the merging black holes 
since the timescale for shrinkage of the binary orbit by gravitational wave radiation becomes 
shorter than the timescale for mass inflow due to viscose stresses in the disk [8]. It is then 
hard to excite transient dissipative processes in the faraway accretion disk. In addition, the 
resulting long times scales of emission correspond more to the afterglow phase than to the 
prompt emission. 

In this paper we take an alternative approach, which follows in spirit the seminal paper of 
Blandford & Znajek [TU]. Blandford & Znajek demonstrated that the rotation of the space- 
time itself can generate a powerful electromagnetic outflow carried by a strongly magnetized 
wind. In case of merging black holes the rotation of the space-time is induced by the orbital 
motion. The properties of the resulting Poynting flux depend only on the black holes masses, 
their separation and the assumed magnetic field, which in turn can be estimated using a 
standard anzats of accretion physics. 

Previously, in Ref. [H] [12] numerical simulations of the electromagnetic signal during the 
last orbit before the merger were performed, showing a production of collimated Poynting 
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flux jets. The jets are produced due to linear motion of the black holes though magnetic 
field via a mechanism that roughly resembles unipolar induction. The classical analogy of 
this mechanism is the motion of planet Io in Jupiter's magnetic field [T3]. The mechanism 
we discuss in the current paper is different, relying on the principles of a Faraday disk, where 
the rotation of space-time plays a role of the conducing surface, drags the magnetic field 
lines and generates electromagnetic outflows [10J. 



II. FARADAY DISK IN ASTROPHYSICS 

Many astrophysical sources are powered by the rotational energy of a central source using 
magnetic field as a "conveyor" belt to transport energy to large distances. Pulsars [H] are 
the prime examples of such process. Qualitatively, the neutron star surfaces permeated by 
the magnetic field works as a Faraday disk or unipolar inductor, generating a Poynting flux 
propagating away from the star. A qualitatively new application of this idea was proposed 
by Blandford and Znajek [ID], who showed that, formally, there is no need for a hard 
conducting disk to launch the Poynting flux: rotation of the space itself can effectively act 
as a Faraday disk. In this paper we apply the ideas of Blandford and Znajek [ID] to the 
space-time produced by orbiting black holes. 

An important difference between the astrophysical applications mentioned above and 
the classical Faraday experiment is that in generating the Poynting flux large electric fields 
develop, which even in the absence of any pre-existing plasma will lead to vacuum breakdown 
due to various radiative effects. Thus, the system forms the plasma out of the vacuum, so 
that in general the problem should be describe by relativistic magnetohydrodynamics. 

The magnetic energy density in the low density regions (from where the outflows are 
generated) is expected to greatly exceeds the plasma energy density, including rest mass. 
Charged particles carry current and generate charge density, that ensures ideal MHD condi- 
tion E • B = 0. In the limit of negligible inertial contribution, the system reaches a so-called 
force-free state, where electromagnetic forces balance itself. We stress that this is not vac- 
uum: inertia-less charge particles ensure ideal MHD conditions. (Formally, in such plasma 
the stress tensor is diagonalizable; in addition, a condition \B\ > \E\ must be satisfied.) 
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III. ELECTROMAGNETIC POWER OF ORBITING MASSES IN EXTERNAL 
MAGNETIC FIELD 



Let us consider two black holes of equal mass M moving on a circular trajectory of 
radius R. Let us assume that there is a magnetized disk outside of the orbit that supports 
the electric currents, which produce a typical magnetic field B within the orbit. We are 
interested in the power produced specifically due to the orbital motion of black holes. If 
black holes are spinning, each will in addition act as a Faraday disk producing the Blandford- 
Znajek Poynting power [TU]. In addition matter accreting onto black holes will produce 
some accretion luminosity both through viscous dissipation in the disk and by launching 
magnetized winds via the Blandford-Payne-type [15] mechanism (as seen in simulations 
Let us give qualitative estimates of the electromagnetic power specifically due to the 
orbital motion of black holes. (Below, for clarity, we omit numerical factors of the order of 
a few). 

The power produced by the Faraday disk can be estimated as 

L EM ~ R 2 B 2 c{Q,r/cf (1) 

where Q is the angular velocity of disk rotation. In the case of orbiting black holes, Q is the 
typical angular velocity of the rotation of the space-time within the black hole orbit. In the 
weak field regime the angular velocity of the rotation of space-time within the orbit can be 
estimated as [e.g. [TBI Eq. 105.20] 

GjM| = (CM) 3/2 
~ c 2 R 3 "' cW 2 1 ' 

where |M| ~ \J GrM 3 / 2 is the total angular momentum of the orbiting black holes (cf. 

Ref. [IT] , problem 13.18). Note, that since outside of the orbit fl oc -R -3 , the corresponding 
electromagnetic power is small. 

Alternatively, we can calculate the rotation frequency of space-time due to orbiting black 
holes of equal masses M on a circular orbit if radius R well before the merger, when the 
motion can be treated in a linearized theory. As a time-averaged approximation, assume 
that two black holes may be represented as a ring of mass 2M rotating with velocity v = 
JGM I (2R), producing a mass current I m = y/GM 3 / 2 /(\/2irR 3 / 2 ). In the linearized theory 
and assuming the time-independent field, the components of the metric tensor satisfy [T8] 

Ago<t> = -167rT O (3) 
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The component then can be found from the known expression for magnetic field produced 
by a current ring [T^l p. 181] 

2\/2 



9o</> 



k 



-(GM/rf 2 



1 = ( ( j- 2 —) K(k 2 ) - ^E(k 2 ' 
Vr 2 + R 2 + 2rRsm6 \\k 2 -lj K J k 2 



4rR sin 9 



(4) 



r 2 + R 2 + 2rRsin6 

where K and E are elliptic integrals. The angular frequency of the rotation of space is then 
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(Note that there is no singularity at r = 0, see Fig. pj). 

In the orbital plane, Q is a slowly varying function of r with the value at r = equal to 
Q = (GM) 3 / 2 /(y/2c 2 R 5 / 2 ), in agreement with our estimate see Fig. [I] 
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FIG. 1. Angular frequency of the rotation of space as a function of radius for a massive ring rotating 
with Keplerian velocity at radius R. Frequency is normalized to (GM) 3 / 2 /(c 2 R 5 / 2 ). Inside the ring 
$7 is nearly constant, falling off oc r~ 3 outside of the ring. Divergence at r = R is an artifact of the 
weak field approximation used to calculate £1. 
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Consider next a Faraday disk of radius R embedded in a highly magnetized plasma in 
the force-free approximation discussed in £jIT] and permeated by a constant magnetic field 
aligned with the disk normal. The disk is rotating with angular velocity Q. The following 
fields expressed in cylindrically collimated r, 0, z then satisfy Maxwell's equation and ideal 
condition E B = [201 121]: 

Qr 

B r — 0, Bj, — E r — — —Bo, B z = B (6) 
for r < R and B z = Bo, E = for r > R. The Poynting flux is then 

L EM = \BlR 2 c{Rn/c)\ (7) 

o 

where VL is the angular velocity of rotation of space ([2]), (and not the angular velocity of the 
Keplerian rotation Qk)- 

We stress that in the model problem considered above (of a rotating disk in a force-free 
plasma) the rotation of field lines at r < R does not break a force balance in the cylindrical 
radial direction r, resulting in a formation of cylindrically collimated outflow. We expected 
that the resulting jets will have a low baryon contamination and thus will be moving at 
relativistic speeds before they start interacting with the surrounding medium. 

Thus, the electromagnetic power due to orbiting black holes acting as a Faraday disk in 
the external magnetic field of strength B can be estimated as 

Lem ~ ^B 2 = 6 x lO^ergs" 1 b\m\ (R/R G y x (8) 
C°K 

where 6 3 = B/10 3 G, m 6 = M/1O 6 M , R G = GM/c 2 = 1.5 x 10 n cmm 6 . (In the numerical 
estimate in Eq. ^ and below the ratio (r/r G ) is treated as a parameter with no dependence 
on mass.) 

(Note, that qualitatively, the power of the Blandford-Znajek process can be similarly 
estimated using Q = ac/r G , 

L BZ ~ a 2 B 2 r 2 G c (9) 

where a is the black hole spin parameters and r G = GM/c 2 is the Schwarzschild radius. At 
the merger moment r ~ GM/c 2 , so that Eq. g agrees with Eq. ^ for a ~ 1, a critically 
rotating black hole.) 

Let us compare the electromagnetic power to the power emitted in gravitational waves 
(GWs). The GW power [TBI §110] 

^ GM 2 n%R* G 4 M 5 

Lgw ? — " (10) 



The ratio of EM to GW luminosities is 

Lem _ B 2 R 4 _ E B 

W ~ GM 2 ~E~ G [ ' 

E B = B 2 R 3 (12) 

E G = GM 2 /R (13) 

where Eg is the energy of magnetic field within the orbit of BHs and Eq is the gravitational 
potential energy. 

The two powers are equal when 

r cq « = 2 x 10 16 cm 63 1/2 mf (14) 



For smaller separations, the GW power dominates. The radius (14) is very large, with a 



merger time ~ c 5 r A /(G 3 M 3 ) larger than the age of the Universe. 



IV. ESTIMATES OF MAGNETIC FIELD 



In the above estimates we used the parametrization of the magnetic field somewhat 
arbitrarily. In this section we estimate the magnetic field within the orbit produced by the 
currents flowing in the accretion disk. As mentioned in the Introduction, two stages of disk 
dynamics may be identified. At early times the loss of energy via gravitational radiation is 
slow enough, so that due to viscous diffusion the inner edge of the disk will be located close 
to the orbital radius, R d « 2R. At later times, for R d < £ d GM/c 2 , £ d « 40 - 100 [8j[22], the 
binary will decouple from the disk, undergoing a merger, while the inner edge of the disk 
remains fixed at R d . The decoupling occurs at time t d ~ £ d GM/c 2 ~ lOyrsmg 3 before the 
merger. 

Several different scalings can be used to estimate the magnetic field in the disk and within 
the orbit. Magnetic field is generated in the disk by the MRI dynamo |23j . which in a thin 
accretion disk produces fields at equipartition with the turbulent motion, B 2 ~ pv 2 urb . In the 
standard a disk the turbulent velocity is a fraction of the local sound speed c s , v tU rb ~ a/«c s , 
while the sound speed c s is a function of the parameters of the central object, mass flow M 
and location in the disk [21] (e.g. in the radiation or matter-dominated parts of the disk). 
It is somewhat beyond the scope of this paper to discuss all the possible combinations of the 
parameters of the accretion disk. Instead, we estimate the upper limits on magnetic field 
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given various parametrizations of the accretion flow. As we will see, even the upper limits 



on the magnetic fields still produce fairly low luminosities, Eq. (19). 

As a first approximation, we assume that the turbulent velocity is of the order of the 
sound speed, v tU rb ~ c s ; m this case the magnetic field is in the equipartition with the 
particle energy-densities, B 2 « pc 2 s . Secondly, we assume thick accretion disk, so that the 



sound speed c s at radius R is of the order of the free fall velocity, c s ~ ^JGM/R. Two 
possible estimates then may be given for the plasma density p. The density can be found 
either from Bondi scaling, pR 2 J 'GM / R w M p ex {GM) 2 j '(? 8jfSX , where p ex and c SfiX are density 
and sound speed in the surrounding medium. This gives 

(GM) 5 ^y/p^ ln 4^ 1/2 -3/2 / R \" 5/4 

b b = ~^rk ~ 6 x 10 4 Gn_ / 1 c se ; 6 — (15) 

where the subscript indicates Bondi scaling and the external number density is scaled as 
p ex ~ m p n = m p (n/0.lcm~ 3 )n_i and the external sound speed as c Sjex ~ 10 7 cms _1 c S)eX) 7. 

Alternatively, the density can be found from the requirement that the accreting matter 
powers an Eddington-type outflow tjmMc 2 = L Edd with some efficiency r\ M « 0.1(^/0.1)^,- 

(GMf/^ _ , . . n _ 1/2 _ 1/2 ( R \ " 5 / 4 



B M = \ ' v " « 3 x 10 5 G m^^l — (16) 

Finally, magnetic field can be estimated assuming that a fraction of the Eddington lumi- 
nosity is carried by magnetic field, B 2 m r]EL-Edd/(cR 2 ), 

Bj .,./«, !xlrt% ^, (*)" (17) 



In all cases the magnetic field within the orbit should be evaluated at the inner edge 
of the accretion disk, located either at Rd ~ R before decoupling or at Rd ~ 100i?c after 
decoupling. Accordingly, we find 

(GM)"/y x f r f 

^ ^4cR fOT GM7J < & 
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The maximal powers, reached at the time of the merger are 

L B ~ 10 37 ergs _1 mgn_ic^ i6 , L M ~ 10 39 ergs" L E ~ 10 38 ergrWte.-i (19) 

Thus, various parametrizations of the magnetic field give a fairly consistent value of the 
electromagnetic luminocity. 

The temporal evolution of the luminosity in each case can be found from the equation 
for the orbital separation at time —t before the merger [T5| §110], 

R = ((GM) 3 (-t)/c 5 ) 1 / 4 (20) 

The total energies emitted after the decoupling are 

E B « v ~"; rux « 10 43 erg 



(GMf Pc 

W^d ^ r ^ 1n 44, 



erf 



(GM) 2 m pV ^ 



E M ~ o ~ 5 x 10 44 erg 

VmC 2 ct t 



{GMfm^E 



E E « ' p ^ — « 5 x 10 4 ^erg (21) 



V. EXPECTED EMISSION 

We foresee two types of electromagnetic signals coming from merging of black holes. 
First, internal dissipation within the jets can produce highly boosted emission similar to 
blazar jets and high energy emission from Gamma Ray Bursts (GRBs). (Blazars are types of 
Active Galactic Nuclei (AGNs) with relativistic jets directed nearly straight at the observer.) 
Blazars dominate the extragalactic sources in the Fermi catalogue [231 125] . This type of 
emission will be highly variable due to relativistic boosting, similar to prompt GRB emission 
[27] . On the other hand, internal jet emission will be highly anisotropic, mostly aberrated 
along the direction of motion of the jets - perpendicular to the binary orbit. Thus, this type 
of emission will be detected only if the line of sight to the binary is nearly aligned with the 
binary normal. 

Secondly, when the jet starts to interact with the circumbinary medium, its kinetic 
power will be dissipated in a termination shock, similar to the case of pulsar wind nebulae 
(magnetically-dominated jets, which lack strong termination shock, may produce similar ob- 
servational signatures [28].) The post-shock flow will be only weakly relativistic, producing 
nearly isotropic emission. 



By analogy with the emission from PWNe, one expects that a broad range of frequencies 
is produced by the Poynting-flux dominated wind. Unfortunately, our understanding of the 
dissipation and accelerations processes in such systems is not sufficient to make quantitative 
predictions. By comparison, e.g. with the Crab nebular, we expect that that the peak of the 
spectral energy density falls into the optical-UV-soft X-ray band, with the total dissipated 
power reaching tens of percent of the Poynting power. 



VI. DISCUSSION 



We have calculated the electromagnetic power expected from the merging black holes. 
The expected Poynting scales not as a square of the orbital frequency, but as a square of the 
rotational frequency of the space-time, Eq. (|2]). The maximum power is reached during the 
final merger orbit, when the orbital frequency and the rotational frequency of the space-time 
are nearly equal, since the resulting black hole has the spin parameter of the order of unity. 
Still, the distinction between the two frequencies is a principal point in our approach. 



Overall, the corresponding peak Poynting powers (19) are fairly low, even though the total 



emitted energies are reasonably high (21). In addition, it is expected that only a fraction of 
the Poynting power is converted into radiation. Still, there is a chance of detection in the 
fortuitous case when the line of sight to the system is nearly aligned with the orbital normal 
and if the outflow reaches relativistic velocities. In this case most of the electromagnetic 
power will beamed along the direction of motion. 

I would like to thank Scott Hughes and Milos Milosavljevic. 
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